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Abstract 

In this paper we construct an effective field theory (EFT) that describes long wavelength grav- 
itational radiation from compact systems. To leading order, this EFT consists of the multipole 
expansion, which we describe in terms of a diffeomorphism invariant point particle Lagrangian. 
The EFT also systematically captures "post-Minkowskian" corrections to the multipole expansion 
due to non-linear terms in general relativity. Specifically, we compute long distance corrections 
from the coupling of the (mass) monopole moment to the quadrupole moment, including up to 
two mass insertions. Along the way, we encounter both logarithmic short distance (UV) and long 
wavelength (IR) divergences. We show that the UV divergences can be (1) absorbed into a renor- 
malization of the multipole moments and (2) resummed via the renormalization group. The IR 
singularities are shown to cancel from properly defined physical observables. As a concrete ex- 
ample of the formalism, we use this EFT to reproduce a number of post-Newtonian corrections 
to the gravitational wave energy flux from non-relativistic binaries, including long distance effects 
up to 3PN (v e ) order. Our results verify that the factorization of scales proposed in the NRGR 
framework of Goldberger and Rothstein is consistent up to order 3PN. 
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I. INTRODUCTION 



Understanding the dynamics of compact binaries within general relativity has become a 
problem of experimental relevance. Gravitational wave detectors such as LIGO/ VIRGO jlj], 
or the planned LISA {2] are capable of probing the evolution of such systems over their entire 
life cycle, from the initial slow inspiral phase to the final ringdown period after the binary 
constituents merge. 

In general, the two-body problem in general relativity is highly non-linear, and therefore 
tractable only by numerical methods. In certain kinematic regimes, however, there is a clear 
separation of scales that opens up the possibility of employing analytical (perturbative) 
techniques to make systematic predictions. One example of this is the post-Newtonian 
(PN) expansion, which is possible whenever the orbital separation r is parametrically larger 
than the typical gravitational radius, of order r s = 2G^m. In that case, one can expand 
the Einstein equations in the velocity parameter v 2 = r s /r <C 1. The PN expansion has 
been thoroughly developed by several collaborations over several decades, reaching a relative 
precision of v 7 (3.5PN) in the computations of two-body motion and of gravitational wave 
observables (energy-momentum and angular momentum fluxes) for non-spinning binaries. 
See j^l for a review and further references. Another limit which is analytically tractable is 
the extreme mass ratio limit, relevant to binaries containing one supermassive black hole. 
This limit can be formulated in terms of black hole perturbation theory, i.e. by finding wave 
solutions to the Regge- Wheeler or Teukolsky equations. See ref. [4]. 

It is interesting to note that analytical control over the two-body dynamics coincides 
with the emergence of a wide hierarchy of scales in the system. Field theoretic problems in 
which there is a separation of scales are known to admit a natural formulation in terms of 
effective field theories (EFTs), and in fact the quantum mechanical counterpart of the (non- 
relativistic) two-body problem in QED and in strong interaction physics (QCD) has been 
successfully addressed using EFT techniques jsj. This motivated ref. js] to recast binary 
dynamics in general relativity as an EFT. 

The EFT approach of jf], Q] attempts to simplify both the computation of observables and 
the physical interpretation of terms in the perturbative expansion by treating each scale that 
arises in the problem one at a time, "integrating out" short distances before moving on to 
the largest scales accessible to experiments. In the two-body problem, the shortest distance 
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scale is the size of the compact binary constituents themselves. Integrating out this scale 
simply consists of writing down a Lagrangian for the center-of mass worldline of each object 
that consists of an infinite hierarchy of terms constructed out of the gravitational field, and 
suppressed by more powers of the object's radius. Given a model for the internal structure 
of the compact object, the coefficients in this Lagrangian are fixed through a matching 
procedure discussed in more detail in ref. jfjj]. We stress, however, that the EFT approach 
is more than just employing a point particle approximation: besides the physical radius 
of the binary constituents, there are several other scales in the problem, and, as discussed 
in jg], properly removing these from the long distance physics necessitates a non-trivial 
decomposition of the gravitational field into modes with support over different kinematical 
regions. 

In recent years, the framework of [6] has developed into a tool used to solve a number 
of problems in gravitational radiation and black hole dynamics within general relativity. In 
particular, the formalism introduced in |(J was extended to account for dissipative effects, for 
example absorption by black hole horizons and tidal heating in neutron stars, in refs. jg, [9|. 
Doing so necessitates the inclusion of additional degrees of freedom living on the black hole 



worldline. Spin interactions were included in the forma 
obtain new results for various spin-spin interactions 11 



ism first in ref. 



101 ]. and used to 



12] within the PN expansion. An 
[[14I which included a treatment 



161 ] at leading order. 



EFT for the extreme mass ratio limit was constructed in [13J, 
of gravitational radiation reaction consistent with standard results 
The calculation of higher order spinless potentials within the EFT formalism, employing 
an extremely convenient parametrization for the metric introduced in refs. 17H19|. was 



initiated in ref 



20), which reproduces known results at order v 4 and paves the way for the 



automatization of terms at v e and beyond. Ref. [21| computed the energy momentum tensor 
and [22j studied higher graviton vertices. Finally, formal applications to the thermodynamic 
phase diagram of compactified black holes (for further references to the literature, see 23j ) 
were developed within effective field theory in 13, |2J, l25j ]. 

In this paper, we consider long wavelength gravitational radiation from compact sources 
within the EFT context. At the linearized level, the relevant framework is of course simply 
the multipole expansion 26J|. In this case, the relevant expansion parameter is a/A 1 
(with a the size of the source), with radiation sourced by an £-pole moment suppressed by 
a relative factor of (a/X) 1 . As has been known for some time, in addition to the multipole 
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expansion, there is a distinct expansion in a parameter r\ = GVm/A <C 1 (with IGjqm < a 
the gravitational radius of the system) that arises from the non-linear nature of general 
relativity. 

While it is straightforward to implement the naive multipole expansion (see for instance 



the textbook |27j). the ^-expansion has a much richer structure [28l43l|. Long wavelength 
(IR) logarithmic singularities arise already at order rj 1 , and at rf there are both IR, as well as 
short distance (UV) logarithmic divergences 1 . In sec. IIII| we set up a systematic expansion 



in powers of 77 in the EFT language of jg]. 

In order to accomplish this, we begin in sec. [TT] by assuming that the short distance 
(multipole) scale a has been integrated out. This results in an effective action consisting 
of a point particle with worldline localized sources (corresponding to the moments of the 
system), coupled to gravity. This effective Lagrangian is manifestly invariant under both 
diffeomorphisms and worldline reparametrizations. In sec. Ill Al we discuss what observables 
can be computed using this Lagrangian. 

Given these preliminaries, we set out to explicitly compute corrections in the post- 
Minkowskian parameter rj, working up to order rf (and therefore order 3PN for non- 
relativistic systems). We find that there are IR divergences at the level of graviton emission 
amplitudes, essentially due to the long range of the Newtonian gravitational potential. How- 
ever, we show that these cancel explicitly from observable quantities, at least to order rf. 
In sec. IIII Dl we consider IR effects at higher orders in perturbation theory. Using a line 



of analysis first developed in |34j, we estimate that the leading IR singularity at order 77™ 
scales like the n-th power of a divergent logarithm. Consequently, we find that summing 
the leading logarithmic IR divergences to all orders results in an overall phase factor in the 
amplitude, therefore canceling from physical observables. This gives partial indication that 
the formalism is free of long wavelength singularities at all orders in the rj expansion. 

In addition to IR logarithms, at order rf there are UV divergent logarithms. Unlike the 
IR effects, the UV singularities are true singularities of the long wavelength EFT, which must 
be dealt with by renormalizing the coupling constants of the theory (i.e., the moments). The 
procedure for doing so is described in sec. IIII CI There we show how the introduction of 



1 The conventional nomenclature for these effects is the "tail" of the gravitational waves [32[ . These have 
been computed in ref. [3o| (at order r/ 1 ), and at order rf in ref. 33 1. 
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renormalized moments renders finite the graviton amplitudes to order rf. (We focus only 
on the quadrupole moment, although the same procedure applies to the higher moments 
as well.). The renormalized moments exhibit subtraction scale dependence governed by a 
non-trivial renormalization group (RG) equation. In principle, this RG equation generates 
and resums the entire series of "leading UV logarithms" at any given order in the expansion 
parameter. Although in practical situations the logarithms do not seem to become large 
enough to necessitate this resummation, the RG scaling is theoretically interesting as it 
gives quantitative constraints on the form of the perturbative expansion. 

The analysis of sec. II III is universal, and applies to long wavelength radiation from any 
localized system. In order to make definite predictions, one must specify the multipole 
moments of the EFT, by matching to a more complete theory that describes the short 
distance gravitational dynamics of the system. Matching is standard in EFTs, but we 
outline how the procedure is carried out in the present context in sec. IIVI For illustration 
purposes, we apply our general framework to post-Newtonian svstems in sec. |VJ There we 
use the results of sec. IIVI together with the NRGR formalism [6], an EFT for gravitating 
non-relativistic systems, to compute the multipole moments needed for the 1PN corrections 
to the energy flux. Together with sec. II I II this allows us to reproduce PN corrections for 
spinless systems at the orders v 2 , v 3 , and v 5 beyond leading order. In addition we compute 
the non-analytic terms v 6 In v at 3PN. Finally, we present our conclusions in sec. IVII 



II. EFT SETUP 



Consider an arbitrary compact source, e.g, a binary system, emitting gravitational radia- 
tion. In this paper we focus on the part of the wave spectrum with wavelength much larger 
than the characteristic size of the source. It is then natural to describe this radiation in terms 
of source multipole moments. Formally, this means that one can systematically decompose 
the motion of the system in terms of a central worldline x^(r) that traces the source's path 
through spacetime, as well as a set of moments that describe the internal dynamics. These 
moments are a set of dynamical degrees of freedom localized on this worldline, labeled by 
their transformation properties under rotations in the rest frame and by transformations 
under parity. 
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In equations, this information is encoded in the Lagrangian 8] 
S = -m f dr- 1 - I dx^L ab ujf{r)+ l -Y^ I dTC^I aba ^(r)V ai • • • V an E ab (x) 

J J n=0 J 

1 °° f 

+ oE / drc^.r b "-(r)V ai ---V an B ab (x). (1) 

Z n=0 J 

Here, dr = ^g fiu dx' 1 dx u is proper time along the central worldline. We have chosen a local 
Lorentz frame at each point x^{r) in such a way that = = dx^/dr is the four- velocity, 
and e^(r), a = 1, 2, 3 comprises a spatial frame whose rotation encodes spin dynamics 2 . Our 
choice of frame therefore satisfies 

r = e ^ Q _ ^ <eV (2 ) 

The components uj ab of the spin connection 3 couple to the total angular momentum L ab (r) = 
-L ba (r) of the compact source. Finally, the SO (3) tensors J a i- a < ; j^u-at > 2), taken 
respectively to be of electric and magnetic parity, and symmetric and traceless with respect 
to the Euclidean metric S ab , define the full set of multipole moments of the compact object. 
For later convenience, we have chosen the normalization constants in the definitions to be 

We will not use any higher moments for the calculations in this paper. The moments J a i - a « ; 
ja,!,. ..at ^£ > 2) serV e as sources of gravitational radiation and therefore couple to the electric 
and magnetic parity curvature tensors, defined in terms of the Weyl tensor as 

E ab = C^ v ^e a a v v el (4) 
B ab = -\e acd ele d v C^ x<T ey, (5) 

and for £ > 2 their covariant gradients, obtained by operating with V a = e^V M , for £ > 2. 
Note that E ab = E ba , 8 ab E ab = 0, and likewise for B ab . Finally, note that in Eq. (JTJ we have 
only kept terms that are linear in the curvature tensor. Terms quadratic and higher, such 
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There is considerable freedom in the choice of local Lorentz frame, and therefore in the choice of variables 
to describe spin. See |1CH12| for further discussions. 

° & is related to the usual spin connection u> AB for the metric vierbein E A=0 '---' 3 (x), g^ v {x) = 
riABE A ll E B y , by the change of basis A a A (r) = e° (t)E» a (x(t)). 



3 UJ 
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as J drE a bE ab , J drB a bB ab , were first introduced in and encode the tidal response of 

the localized system to external gravitational fields. 

By construction, Eq. ([I]) provides a fully diffeomorphism invariant description of a com- 
pact source, in the limit of long wavelength radiation. Together with the Einstein term 
Seh — ~^ m Pi J d^Xy/gR, with m 2 Pl = 1/(32ttGn), Eq. (JTJ can be used as a starting point 
for a systematic long wavelength expansion, both in the multipole expansion parameter 



a/ A <C 1 and the "post-Minkowskian" (borrowing the terminology of 29]) parameter 
i] = GnM/X <C 1 that controls non-linear corrections. 

In order to compute these corrections systematically, one must specify (as in any other 
EFT) a power counting scheme for bookkeeping relative sizes of terms. Since we are inter- 
ested in radiation of wavelength characterized by a scale A, we assume that derivatives acting 
on the gravitational field scale as ~ A -1 . Hence, we take ~ A, and by demanding that 
the kinetic term for h^ u /m P i = g^ u — r]^ be leading order, h^ u ~ A -1 . The £-th multipole 
moment is expected to scale as ma e , and thus multipole couplings in Eq. (CQ) are suppressed 4 
by powers of a/A <C 1. With these rules, it is possible to determine the scaling in powers of 77 
or a/ A of any term in the action. For instance at a fixed order in a/A, the post-Minkowskian 
terms are generated by Feynman graphs with increasing numbers of insertions of the mass 
monopole term — m J dr ~ mX and n-graviton vertices, which scale as (mpi\) 2 ~ n . 

The rest of this paper is devoted to explicitly computing corrections to low frequency 
gravitational wave emission, both in powers of rj and in powers of the multipole expansion 
parameter. In the next section we give a general discussion of how such observables are 
computed in practice, given a prescribed set of moments / ai --- a ^(r), J ai - ae (r). In sec. Ill II we 
use this Lagrangian to compute corrections to gravitational radiation suppressed by powers 
of rj, working up to order rj 2 . Finally in sec. IIVI we discuss how to determine the moments 
jai...a^ ja x ...a t £ Qr a i oca ij zec i gravitational system. We apply these results to post-Newtonian 
systems in sec. [V] 



4 In practice, for instance in applications to PN systems, the £-th multipole scales not as a definite power, but 
rather as a series ma 1 [l + 0{a/X) + O(afX) 2 + •■•]. Strictly speaking, to have manifest power counting 
one would have to introduce a separate coupling in Eq. (p} for each term in this series. This would be 
extremely awkward to write out explicitly, so we will stick with the standard conventions in the general 
relativity literature of classifying moments by their SO(3) transformation properties rather than by their 
scaling with respect to the expansion parameter. 
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A. Calculating observables 



As in our previous work, we find it convenient to set up diagrammatic rules for the 
computation of observables associated with the emission of gravitational radiation (e.g., 
energy, momentum, and angular momentum flux). In this paper, we will restrict ourselves 
to time averaged observables (a field theoretic formulation of instantaneous quantities is 
given in [14J.) A natural way to set up the diagrammatics is to formulate all observables in 
terms of the matrix element for the emission of a single graviton from the compact source 
described by Eq. (0Q), 



iA h {k) = . (6) 




Here Ah{k) denotes the probability amplitude to emit a graviton of momentum k and definite 
helicity h = ±2 (as measured in a nearly Lorentz frame infinitely far from the source). This 
is given by Feynman diagrams with one on-shell external graviton (internal vertices and 
propagators are obtained from the source terms in Eq. (CQ) plus the Einstein-Hilbert action 
Seh)- 

In terms of Ah(k) one can compute a polarized graviton emission rate as 

1 r/ 3 k 

dlVk) = - - " „ , |A(k)| 2 , (7) 
v ; T(27r) 3 2|k|' v n ' v ; 

where T — > oo represents total integration in the detector, and drops out of time averaged 

quantities. Moments of the differential rate <fT/j(k) give rise to physical observables. For 

example, the (polarized) total rate of radiated linear four-momentum is given by 



h=±2 



k»dr h (k), (8) 



where = (|k|,k) is the four-momentum of the emitted graviton 5 . It is also possible 
to compute the rate of angular momentum loss by the system. In terms of the helicity 
amplitudes Ah(k), it is 

j = X) / hndr h(k), (9) 
h J 

where n = k/|k| is the direction of the emitted graviton and we sum over the physical 
helicities h = ±2. 



5 Strictly speaking, there should be a cut in the integration over graviton frequency at a value w» ~ 1/a 
where the multipole expansion begins to breaks down. 
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As an example, ignoring post-Minkowskian corrections, Eq. ([I]) gives for the one-graviton 
matrix element, 



jijk 



iA h (k) 



(10) 



where we work with physical on-shell graviton degrees of freedom with polarization tensor 
e,j(k, h) satisfying the usual conditions k l ejj(k, h) = e^k, h) = 0. Here we have assumed for 
simplicity that the source is at rest, and take the frame {e°, e a } to align with the Lorentz 
frame at infinity (in which case the rotation of the system is reflected in the time dependence 
of the moments). Then using the standard result for the sum over graviton polarizations 



I 

2j%(k,/i)e* s (k,/i) = - 



Si r Sj S ~\~ S{ S Sj r SijS rs -\- ^2 (^2jk r k s -\- $ rs k^kj) 



(^irkjk^ + 5j s kjk r -|- 5j r kjk s -|- 5j S kjk r ) + ^kjkjkj.k s 



one finds for the total radiated power 



p c 



G 



N 



TFT 



dk 



k e o 1 fi ok 8 9 

- \r^k)\ 2 + —k 6 \r{k)\ 2 + — \i« k {h)\ 2 + ■ 



Transforming to the time domain, one then gets the standard result [2 



,(H) 



(12) 



\dt 3 ' 



+ 



16G 



N 



45 



dt 3 



J ij {t) 



+ 189\V^ 4 



+ 



(13) 



where the time average is (• • • ) = ^ J Q T • • • , as T — > oo. 



III. POST-MINKOWSKIAN CORRECTIONS 



At the linearized level, long wavelength radiation from a localized source is described 
by Eq. ( flOl) . We now show how to compute corrections in the EFT arising from graviton 
self-interactions. Not only are these corrections relevant for gravitational wave phenomenol- 
ogy, but they also introduce logarithmic singularities at both long and short distances. In 
this section, we discuss the physical origin of these divergences and explain how they are 
systematically removed from observable quantities. 
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A. Corrections at order r/ 1 

Long distance (infrared) logarithmic singularities appear already at order rj 1 , from the 
graph in Fig. [TJ This graph represents the interaction of outgoing gravitational radiation 
(sourced, for example, by the system's electric quadrupole) with the Newtonian potential 
generated by the total mass of the system. Using the Feynman rules derived in [6|], we find 
that after projecting onto transverse traceless parts and reducing tensor integrals to scalars 



via standard methods (see e.g., [351]). the graph is proportional to 

iA v i = iA v o (G N mk 2 ) x I(k), (14) 

where 



iA, = ^el^h)Pi(\k\), (15) 



is the leading order quadrupole radiation amplitude, and the function I(k) is a linear com- 
bination of integrals of the form (with d denoting the dimension of spacetime) 

1 V f d d - l a f 1 \ 1-n 1 

W q 1 M (16) 



k 2 J J (2vr) d - 1 WJ k 2 - (k + q) 2 + k' 
with e + . 

The general structure of /(k) is not difficult to understand. Focus on the term with 
n = 0. Then the factor of 1/q 2 is the Fourier transform of the Newtonian potential of the 
source, while the factor 1/k 2 — (k + q) 2 is the propagator for a graviton of energy |k| emitted 
from the quadrupole vertex. Note that in d = 4 dimensions, the integrand for n = has the 
long distance (q — > 0) behavior 

* ,<5 1 1 (17) 



(2tt) 3 q 2 k • q 
10 



and is therefore infrared logarithmically divergent. This IR singularity is physically due to 
the interaction between the nearly on-shell emitted graviton and the long range 1/r potential 
of the source over an infinitely large amount of time. It is entirely analogous to the divergent 
phase factor that appears in quantum mechanical Coulomb scattering. 

In order to regularize this divergence, we use dimensional regularization, keeping the 
spacetime dimension d arbitrary in the calculations and analytically continuing to a neigh- 
borhood of d = 4 in the complex d-plane. With the aid of the table of integrals in appendix |A~1 
we find for the 0{r] 1 ) amplitude 



iArf- — iA v o x {iG^m\\s\) 



(k 2 + ze) 

7T/X 2 



o7e 



(d-4)/2 



X 



2 11 / j t . /vr 2 203 \ 

^4-y + (rf - 4) U + 72-) + 



(18) 



where terms of order (d — 4) 2 and higher have been dropped. 

The scale \x is an arbitrary subtraction scale put in by hand to make the expression 
dimensionally correct as d — >• 4. Formally, it arises by introducing a "renormalized" Newton's 
constant, 

- 2m* 2 J d 4 x^R = -2Z G V)™ V~ 4 / d*Xy/gR, (19) 

or Gfj = Z G {ji)ji A ~ d GM ) with renormalization factor Z G = 1 + 0{h). The \i independence 
of the bare theory then implies a (classically trivial) renormalization group equation for the 
renormalized Newton's constant, 



d 

djj, 



(d — A)G 



N ■ 



(20) 



which will play a role later in Sec. IIII CI From now on, only the renormalized Gm will appear 
in our calculations. 

At order rj 1 , the /i dependence explicitly drops out of physical quantities, however. This 
follows since at this order, the squared modulus of the emission amplitude is 



A 



and from Eq. ( fl8l) . 

An^ 



Re 



A^ 



(GVm|k|) 



A r p 



(d-4)/2 



1 + 2Re4^ + 0( V 2 

AqO 



(21) 



sin ii 



d-4 
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d-4 



+ (d - 4) 



vr 2 203 \ 



Hence as d — > 4 along the real axis, 



A 



A q o 



1 + 27rGArm|k| 



(22) 
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P 3 m m 



P 3 m m 



P 3 m m 




(a) 



(c) 



FIG. 2: Second order post-Minkowskian corrections. 



Although we have only presented this result for quadrupole emission, it is in fact universal, 
and applies to radiation from any multipole moment 6 . Note that, upon squaring, Eq. (|18|) 
also produces infrared single and double poles in d — 4 at order rj 2 . It is an important check 
of the formalism that these poles cancel from physical quantities, as we now show. 

B. Effects at order rj 2 

The order rf corrections to quadrupole emission are shown in Fig. [2j In addition to 
infrared divergences similar in nature to those encountered at order rj 1 , a new feature at 
this order is the appearance of logarithmic short distance (UV) singularities, whose physical 
origin and resolution will be discussed below. 

Operationally, the presence of UV divergences can be seen by examining the structure of 
the graphs in Fig. |2j Take for example Fig. EJ^b), which is proportional to the integral 



where the tensor structure of has been suppressed for clarity. By dimensional analysis, the 
integral over p scales as (l/q 2 )( 5 ~ d )/ 2 , and in coordinate space reflects the interaction of the 
emitted graviton with the [G^m/r) 2 relativistic potential (in d = 4) of the source. (This 
is particularly clear in Fig. EJ^a), which contains a subgraph with two m insertions and one 

6 This follows from the fact that this term is determined by the coefficient of a logarithmic IR singularity. 
We have, in fact, explicitly checked that Eq. (|22l) also holds for magnetic quadrupole and electric octupolc 




(23) 



radiation. 
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internal three-graviton vertex. This is precisely the Feynman diagram corresponding to the 
1/r 2 potential.). From the scaling of the p integral in Eq. ([23]) . we deduce that the short 
distance (q — > oo) behavior of Fig. [2^b) is, in d — A, 

d 3 q 1 1 f rf 3 q 1 



(27r) 3 |q|k 2 -(k + q) 



(2 



7T 



|q| 



(24) 



Fig. MJo) is therefore logarithmically ultraviolet divergent. On the other hand, the integrand 
behaves as l/(|q|q • k) in the infrared limit q — > 0, so Fig. E](b) is infrared finite by power 
counting. A similar estimate reveals identical UV and IR behavior for the graph in Fig. [S^a). 

Performing a chain of reductions to scalar integrals, projecting onto transverse traceless 
polarizations, and employing a set of standard Feynman integrals collected in appendix [Aj 
we find 

(k 2 + ie) i (d " 4) 



A, 



(a) 



and 



A. r p 



89 

^(G N m\k\y 



-(G N m\k\f 



7T fX z 

(k 2 + ze) 

Tlfl 2 



(d-4) 



X 



1 14989 
d-A ~ 18690 

1 157' 
d-A ~ To" 



(25) 



(26) 



In these expressions, poles in d — 4 denote the purely ultraviolet divergences discussed above. 

Evaluation of the graph in Fig. Wis) requires more care, as it contains both UV and IR sin- 
gularities. Suppressing the tensor structure from the three-graviton vertices, the propagator 
structure of Fig. Etc) is 

d d-i q d d-i p ill i 



4^(k) 



(27r) d - 1 (27r) d - 1 q 2 k 2 - (k + q) 2 + ie p 2 k 2 - (k + p + q) 2 + ie ' 



(27) 



which represents the amplitude for the emitted graviton to interact twice with the source's 
Newton potential. Infrared divergences come from the region of integration q, p — > 0. In 
this limit, the q, p integrals factorize, and thus 



d 3 q 1 1 



(28) 



(2tt) 3 q 2 k ■ q 

In d dimensions, this produces an infrared double pole l/(d — 4) 2 which, as we will see 
explicitly below, exactly cancels the infrared double pole from the square of the amplitude 
Arji computed above. 

In order to power count the UV divergences of diagram Fig. [2](c) we need to be more 
careful about the dependence of the numerator on the momenta. The diagram contains 
two insertions of graviton self-interaction vertices, each scaling as two derivatives. This 
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introduces in total four powers of momenta in the numerator. Of all the terms in the 
numerator, we focus on a term in the amplitude proportional to k 2 q 2 . For this component, 
the behavior of the integral in the UV region is now 

d d_1 q 1 f d d ^p 1 1 



(2n) d - 1 q 2 



(29) 



(2tt)' j - 1 p 2 (p + q) 2 ' 

Given the |q| rf_5 scaling of the p integral, for d — 4 the q, p — > oo behavior is UV logarith- 
mically divergent, 

^ k > - / OT- < 3 °> 

The exact calculation of Fig. W(c) (or "mother of all tails") is somewhat involved. We 
have found it useful to apply the Mellin-Barnes method for calculating multi-loop Feynman 
integrals, as reviewed for instance in Ref. 35] . Leaving the details of the calculation for 
appendix [A] our result for Fig. |2f c) is 

(d-4) 



A {c) 

-f- = (G N m\k\f 



(k 2 + ie) 

7lfi 2 



olE 



X 



+ 



109 1 



7vr 2 1364777 



(d-4) 2 315rf-4 
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132300 



(31) 



where, as discussed above, the double pole in d — 4 indicates an IR divergence. On the 
other hand, the single pole contains both UV and IR logarithmic singularities. It is useful 
to explicitly isolate the UV pole. This is done in appendix El with the result 
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uv 
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so we can write 
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The total amplitude at order t] 2 is then 
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To compute physical quantities at order i] 2 , one needs the full modulus squared, including 
terms from the square of the amplitude A v i calculated in the previous section. The full set 
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of terms at 0(rf) is 



and from Eq. (jig) . Eq. (l3"3]h 
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Comparing the two expressions, we see the explicit cancellation of infrared divergences. Our 
final result is 
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1 + 2vr (G N m\k\) + (GArm|k| 
+ 0( V 3 ). 
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3 44100 _ 

Although the infrared divergences have disappeared to 0(rj 2 ), the ultraviolet divergence 
remains. In the PN context, the presence of such logarithmic singularities was first pointed 
out in ref. {^I. See also j^. This UV divergence represents a true short distance singularity 
of the EFT, and requires renormalization of the theory as we now discuss. 



C. Renormalization 

It is straightforward to identify a suitable counterterm that cancels the UV pole in 
Eq. (|34|) . From the form of the leading order quadrupole amplitude, we can interpret the 
frequency domain quadrupole moment as a (frequency dependent) bare coupling constant. 
Then the explicit regulator dependence of the amplitude cancels that of the bare coupling, 
leaving behind a finite remainder. 

To make this more explicit, we begin by introducing a renormalized quadrupole moment, 

I*( U ) = Z{u,ii)I*(LJ,ii) (39) 

where we denote the wave frequency u> = |k|. Note that in order to keep the bare moment 
1$ scale independent, we must introduce \x dependence in the renormalized moment 1$ If 
the renormalization constant Z(a;,/i) is adjusted as 

1 



i nv 

Z MS (u^) = l + — (G N muj) 2 x 



(d - 4)uv 



+ 1e — hi47r 



(40) 
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then Eq. ( 138|) becomes 

2 



A 



AR 



1 + 2tt (Gjvmo;) + (GnitluY 



214, k 2 4vr 2 634913 



In — - H h 

105 4/i 2 3 44100 



+ 0(77 3 ), (41) 



where now A^ is the leading order quadrupole amplitude, written in terms of the renormal- 
ized quadrupole I^(uj). Note that, being a physical quantity, \A\ 2 should be independent 
of the arbitrary renormalization scale fi. Indeed from the the fi independence of the bare 
quadrupole Z^(w), together with Eq. (1201) . Eq. (1401 . one finds a renormalization group (RG) 
equation 



//— = - — (G N mu) 2 I*(u,n), (42) 



-f§(GivrM 2 

Iij(uj,iJ,o)- (43) 



d rB , , 214 

and thus from Eq. (138]) . the amplitude is /i independent, at least to order rf. Alternatively, 
Eq. (142]) could have also been derived less formally, by forcing the result of Eq. ( 138]) to be 
scale independent, and introducing fi dependence of the renormalized moment in order to 
compensate for the explicit scale dependence of the logarithm. 

The RG equation for the quadrupole moment has the solution (dropping the superscript 
R from now on) 

J -(w,//)= — 

Because of RG invariance of physical observables, any value of /i can be used in calculations. 
In practice, it is useful to choose fi to minimize the logarithms in \A\ 2 , so one can take /i ~ u, 
where u is the typical graviton frequency in the problem. The scale no should be taken to be 
of order the typical short distance scale, e.g. «o ~ a -1 , where a is the size of the system (the 
parameter that controls the multipole expansion). At the level of the multipole EFT, /i is a 
free parameter, to be fit to data, and dependent on the details of the underlying gravitating 
system. However, if the physics is known at the scale /io ~ where the multipole EFT 
begins to break down, it is possible to do an explicit matching calculation to fix the precise 
value of Un. This is a standard procedure in effective field theories, see for example the 
reviews [37j . We will come back to this issue in a separate paper. 

In typical applications, the renormalization group is used to sum up the series of large 
"leading logarithms" that arise whenever log /xo/a* becomes large enough to compensate for 
the smallness of the expansion parameter, thereby improving naive perturbation theory. 
Unfortunately, in gravitational wave physics the logarithms cannot become large. This is 
because in the logarithmic terms r] 2 hifj, /fi ~ rf In a/A with r\ ~ r s /\, the extent of the 
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gravitational wave source is bounded to be a > r s and therefore the logarithm can never 
compensate for the smallness of rj. For example, in the PN regime, one expects aWa 2 ~ v i 
while the post-Minkowski expansion parameter becomes i] ~ v 3 , so that rf In Ho/fi ~ v®\nv 
is still small. In this case, terms higher order than [G^mu) 2 In a*o/a* i n Eq. (|43|) are no 
larger than uncomputed PN corrections, and need not be resummed. Nevertheless, the re- 
summation of UV logarithms from the RG is systematic and can be used to obtain resummed 
gravitational wave observables in a factorized form, including instantaneous observables such 
as the waveform. 

While summing the full series of logs in Eq. ( T4"3"j) is typically not necessary for phenomenol- 
ogy, the RG equation for the quadrupole does contain information about the dynamics, as 
it constrains the pattern of logarithms can appear in the amplitude squared at higher orders 
in the calculation of quadrupole radiation from compact systems, 



Ao(w,/i ) 



1 428,^ , 2l u , 91592.^ 4 2 w 
= 1 - —z {G N mw In — + — — {G N mw In 2 — 

leading log 105 /*> 11025 fl 

39201376.^ , 6l 3 u /AA . 

"ImtW^ ln V'-'- (44) 

This set of terms is universal, independent of the short distance dynamics. It is a prediction 
of the RG that this series of terms will appear in long wavelength gravitational radiation 
from any system. To our knowledge, this series of terms, with the precise coefficients 
given above, has not been given previously in the literature. It should be possible to test 
this prediction in the extreme mass limit, by expanding out the solutions of the Regge- 
Wheeler equation to the required order. While Eq. ( l4"3j) does not seem to have immediate 
phenomenological applications, it may be possible to use this equation, for instance, to 
construct phenomenological templates for gravitational wave emission, along the lines of 
refs. [38]. 

We have referred to the series of logarithms in Eq. ( ]44"1) as the "leading logarithms" . This 
means that this series captures the set of UV divergences of the form rf n ln n A, n > 1, where 
A — > oo is a UV cutoff on integration momenta. By power counting arguments similar to 
those applied in the previous section, one can show that there are new UV divergences, 
proportional to In A, at every even order in the rj expansion. These divergences can again be 
rendered finite by a renormalization of the multipole moments, with the consequence that 
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the RG equation Eq. (j42p has the following structure to all orders in rj, 

( 45 ) 

where /?2 = —214/105, and /32n>4 is obtained by computing UV poles in diagrams with 2n 
mass insertions. Formally keeping only the contribution due to 02 then sums the leading 
logs, keeping /?4 fixes the coefficients of all next-to-leading logs of the form rj 2n ln n_1 A for 
all n > 2, and so on. Note that this discussion neglects UV divergences in Feynman graphs 
with insertions of spin couplings. These may introduce additional UV poles at odd orders 
in n. 



2n 



71=1 



D. IR structure at higher orders 

We have only shown the explicit cancellation of IR divergence up to order rj 2 , but clearly 
it should persist at higher orders. A relatively simple class of such divergences is presented 
in Fig. [3j These graphs contain the most singular IR divergence at each order in rj. We have 
explicitly computed the first two of these graphs in the previous section and shown them to 
have l/(d — 4)j_r and l/(d — 4)^ R divergences. It is not difficult to show that this pattern 
persists at higher orders. 

Consider a graph with similar topology to those in Fig. [31 with n insertions of the mass 
monopole. Such a graph contains integrations over n spatial momenta qi=i...., n - By examin- 
ing the region of momentum where all q^ — > 0, one finds the behavior 



iAn iA v o x — 
nl 



-32n(G N mk 2 



d d ~ L q 1 1 
(27r) a! - 1 q^2k- q_ 



(46) 
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Then, using the result 



d^q 1 



(27r) d - 1 q 2 2k • q - ie 
one finds that the leading IR poles sum into 



16vr|k| (d-A) IR 



(47) 



A n — > A n o exp 



n=0 



(d-4) 



(48) 



The IR divergences have summed into a harmless overall phase factor, and therefore cancel 
from physical quantities. The summation of IR divergences due to Coulomb-type phases in 
soft graviton emission processes was first discussed by Weinberg j^ . 

Associated with the leading IR divergences are additional finite parts, which do contribute 
to |*4.| 2 . These contributions come with factors of it at every order in rj, and tend to be 
numerically enhanced. The precise form of these corrections can be obtained by solving the 
wave equation for the propagation of a graviton in the presence of the 1/r gravitational field 
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of the source, see refs. 
1*4.,/) | 2 gets scaled by the Sommerfeld factor 



40j. Their result is that the leading order amplitude squared 



1 — exp (— AttGntuw) 



(49) 



411 ] in the construction of phenomenological wave- 



This result has also been employed in 
forms. 

Taken together, the resummation of the leading UV logarithms of Eq. ( 1441) and the 
resummation of the leading IR effects via Eq. fj4"9l predict the existence of a term at order 
rf that has the form 

D — T^{G N muY\n— . (50) 



AqO 



8567T 3 U) 

105 fio 



This prediction agrees with results for extreme mass ratio post-Newtonian systems, calcu- 
lated by the methods reviewed in 4J]. 



IV. THE MULTIPOLE EXPANSION 



It remains to match onto the multipole EFT of Eq. (CQ), that is, to compute the specific 
form of the moments given the details of the short distance physics. In order to match, one 
must compute graviton emission amplitudes in the full theory (containing details about short 
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distance structure), expanded in the long wavelength limit, and compare them to Eq. ([I]). 
One then adjusts the moments in Eq. ([1]) in order to obtain agreement between the two 
calculations. 

For classical processes, it is often sufficient to consider single graviton emission ampli- 
tudes. These are generated in the full theory by the interaction 

T[h] = -— L- / d A xT^(x)h^ (51) 
zmpi J 

where the tensor T^ l/ (x) should be regarded as a Minkowski spacetime tensor, measured 
relative to a Lorentz frame asymptotically far from the source. T^ v (x) is defined in terms 
of the off-shell amplitude for single graviton emission 

A- = -7T^— I d*xT^{xy k - x , (52) 
zmpi J 

and therefore satisfies the conservation law d p T^ l '(x) = on account of the Ward identity for 
graviton amplitudes. The tensor T^ u (x) can therefore be regarded as the energy-momentum 
"pseudo-tensor" for gravity plus matter that appears in all approaches to gravitational 
radiation. 

The multipole expansion is generated by taking the limit k — > 0, or in coordinate space, 
by expanding Eq. ( 15"TT) as 



oo 1 

M*) = E ■ ■ ■ ^ d h ■ ■ ■ di n M*°> °)> ( 53 ) 

n=0 

where the point x = is taken to be the center of mass, defined such that 

X* = / d 3 xT 00 x ! = 0. (54) 



This expansion can be truncated at finite order whenever the condition x - V ~ a/A < 1 
holds. Plugging into Eq. ( )5T|) and breaking up the moments 

c^xT^x* 1 • • • xS (55) 

into irreducible representations of the rotation group gives the result Eq. flTJ, at linear order 
in the radiation field h^ u . The non-linear terms, usually not needed for classical applications, 
then follow by covariantizing the linear result. 

The decomposition of moments of T^ u into irreducible representations is standard in the 



general relativity literature, see for example 
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27J. Here we recall some of the manipulations 



at low orders in the expansion. To zeroth order in derivatives, the multipole expansion just 
gives 

r[7j] = -_ L_ / dx°m/i o(x ,0), (56) 



2m 



PI 



where m = J d 3 xT 00 is the total energy of the isolated system. The relation d^T^ u = 
implies that this is a conserved quantity, m — 0. Eq. (156]) . of course reproduces the mass 
monopole term — m J dr in Eq. ([1]), to linear order in the fields. We have not included in 
Eq. (156]) the contribution from 



dx c 



d 3 xT 0i ( n x 



(57) 



since by conservation of T Ml/ , 



rf 3 xT 0i (x°,x) = / ^xT uu (/,x)x ! = mX' = 



3 rfOO,' 



(58) 

in the center of mass frame. Even if this term were non-zero, for on-shell radiation do ~ <9j, 
so this coupling to h 0i would not appear until first order in derivatives. Likewise the coupling 



dx° 



d 3 xr j (x°,x) 



hij(x°,0), 



does not appear until second order in derivatives, given the relation 



d 3 xT iJ '(x°,x) = -dl d 3 xT 00 (x ,x)xV. 



(59) 



(60) 



which follows from energy-momentum conservation. 

Including all terms that survive at first order in derivatives, one finds 



The linearized spin connection is 



to 



2m 



•i / dx*L ij uj^. 



(d%-^), 



pi 



(61) 



(62) 



and the angular momentum is given L* J = — J rf 3 x(T *x J — T 0j x*). This agrees with the 
terms — | J dx^La^uj 1 ^ in Eq. ([T]), again to linear order in the fields. Note that by d^T^ = 0, 
the angular momentum is conserved Lij = 0. 

Multipoles which source radiation start appearing at second order in derivatives. At this 
point, it is convenient to split up the different multipoles into radiative terms and non- 
radiative terms. The radiative terms couple to hij and therefore source physical radiation. 
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The non-radiative terms are couplings either to constants of the motion, as in Eq. (1561) . 
Eq. flfJTj) . to quantities that can be made to vanish by a choice of inertial frame, or couplings 
to the Ricci tensor and its covariant derivatives. On account of the vacuum field equations, 
Rnv — 0, the latter can be made to vanish by field redefinitions. 

From now on we will omit all non-radiative terms. To work out the radiative couplings, 
it is sufficient to set Jiqq = Jiqi = and retain only the physical radiation field hij. As 
discussed above, there is one term at second order in derivatives of the form (TF stands for 
the traceless part of the tensor) 

TF 



T[h} 2 = ~jdx° J d 3 xT 00 xV 



Eij, (63) 



with the linearized electric parity Weyl tensor for on-shell 7 h^, 

There are two additional contributions at second order in spacetime derivatives. One 
arises from the expansion of hij to first order in spatial gradients, and is proportional to the 
moment J d 3 xT y x fe , which decomposes into 5*0(3) irreducible representations as 

(2©0)<g>l~3 jB ©2 M ©l©l. (65) 

The representations 2, 3 correspond to magnetic quadrupole and electric octupole moments, 
while the two 1 representations are non-radiative, and vanish by the use of the equations 
of motion. Of the radiative terms, the 2 contributes at two-derivative order. Carrying; out 



standard manipulations which will not be repeated here (see, e.g., the textbook }27| ) one 
finds a contribution 



d 3 xT ij x. k d k h 4 



'k'Hj 

where we have defined the leading order magnetic quadrupole 



3 



~J ij Bij, (66) 



Jij = -lj rf3 x (^ikl [ T °k x j x lfF + e j M pOk^TF^ (6?) 



The magnetic parity Weyl tensor is 



I i 

■^ij "^^imnC'Qjmn ~Z ^imn (^O^n^jm 9jO m hQ n ^ . (68) 
Z ZTTlpi 



7 Had we kept terms with hoo and hoi we would have found the same result, but with the full expression 
Eij — > (dodihoj + dodjfiQi — didjhoo — d^hij) / (2mpi), as expected from gauge invariance of the action. 
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The remaining contribution at two derivatives is generated by the term in the expansion of 
hij containing two spatial gradients. This couples to J d 3 xT y x fe x J , whose decomposition is 

(2 © 0) g> (1 g> 1) 8 ~ 4 S © 3 A / © 2 E © 2 © 2 © 1 © © 0. (69) 

It turns out that the last five moments couple to the Ricci curvature, and are therefore 
irrelevant for radiation. However, the term 



-\T ik yjx k - iT^x'x*] TF <9^% (70) 



has parts which contribute to the two-derivative action. One can simplify Eq. ( ITU]) by 
applying d^T^ = 0, together with integration by parts, with the result 

~ J d 3 xrW^M*>)] = \jd^ {r kk - |r°* x fc + ^r 00 x 2 ) [xvf^, 

(71) 

This contains the second derivative term we are after, together with higher order three- and 
four- derivative terms. Putting everything together, we find at second order in derivatives 

Y[h} 2 = lfdx° (FEv - P^B^ , (72) 

with 

I ij = j rf 3 x(T 00 + T kk ) [x*x»] TF + 0(d l ) } (73) 

and J u given in Eq. (|67l) . 

The pattern is similar at higher orders. At third order in the derivative expansion, there 
are contributions from the 3^ in Eq. (1651) . from 3m in Eq. ( 16 9 j) . and from the coupling to 
j d 3 xT y x rx m , which induces order d 1 corrections to the magnetic quadrupole and one of 
the leading parts of the electric octupole moment. The terms in the quadrupole moment 
at order d n arise from the expansion of hij in spatial gradients at orders n + 2, n + 1 and 
n. More generally, the set of terms in an electric parity £-pole moment at order d n receive 
contributions from the expansion of hij in spatial gradients at orders n + £, n + £ — 1 and 
n + £ — 2. Magnetic parity £-poles at order d n receive contributions from the expansion of 
hij in spatial gradients at orders n + l+1, n + £ and n + £ — 1. 
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Summarizing our results at low orders in the expansion, we find for the first three moments 
which can radiate, 

nh] = \j dx° (piEi, - p^B^ + h l3k d k E l3 + ■ ■ ■ V (74) 

with 

F = J d 3 x (t 00 + T kk - ^T°V + ^T 00 x 2 ^) [xV] TF + 0(d 2 ), (75) 

J* 3 '* = J rf 3 x(T 00 + T") [xVxf^ + O^ 1 ), (76) 

jij = ~J ^xfe^^Vx' + e^r"^) + £>(d 1 ). (77) 

Given a model for the gravitating source (i.e., a prescription for computing the off-shell 
emission amplitude T^ u ), these formulas and their generalization to higher orders in deriva- 
tives can be used to compute long wavelength radiation observables. The results presented 
here are sufficient for computing post-Newtonian radiation to order v 2 beyond the leading 
order quadrupole formula, as we discuss in the next section. 



V. APPLICATIONS TO POST-NEWTONIAN SYSTEMS 

In order to apply the general formalism developed above to post-Newtonian systems, 
one must compute the amplitude T^ v (x) in terms of the orbital degrees of freedom of the 
non-relativistic compact system. We will only consider spinless systems in this paper and 
leave the issue of spin for future work [421 ]. 

In momentum space, T^ v (k) can be calculated as a sum over Feynman diagrams with a 
single external (off-shell) graviton of momentum h** '. By definition, external propagators are 
stripped off. For PN systems it is convenient to perform this calculation using NRGR jfj], 
which is tailored for non-relativistic gravitating bound states. In this EFT, the gravitational 
degrees of freedom are decomposed into potential modes H^ v with support over the size r 
of the bound state, and on-shell radiation modes h^ u . The NRGR radiation mode exactly 
coincides with the graviton field appearing in the multipole EFT written down in Eq. (CD). 
The potential mode, being a short distance field, must be integrated out. In diagrammatic 
language, this means simply that T /iI/ (fc) is computed by considering only graphs that contain 
internal potential graviton propagators, with a single external radiation line of momentum 



24 



FIG. 4: Feynman diagram yielding the leading terms for T 00 and T°\ 
/c M . The relevant Feynman rules are given in 

aa. 

In particular the propagator for the 

potential graviton is 

(H^x^H^O)) = -^(2n) 3 5( q + k)5(x°)P^ (78) 



where in the gauge employed in P^afi = \ [Vi^V^P + VvPVva ~ -^V^Vafi] ■ 

To illustrate how this works in practice, we now compute T fJ,u (k) to next-to-leading order 
in the velocity expansion. We stress, however, that the procedure is well defined to all 
orders in the expansion parameter v. To be definite, we consider an ensemble of gravitating 
point particles with masses m a and positions x (i), with a labeling the particle species. 
Introducing the partial Fourier transform T 7 ^ (a; , k) = J <i 3 xe _ * kx T MI/ (a; , x), we can read 
off the moments from the Taylor expansion about k = 0, 

T"V, k) = ^f- ( [ d 3 x.T^(x°, x) x 41 • ■ ■ x^ k n • ■ ■ k ln . (79) 

n=0 ' ^ ' 

By construction, in NRGR, each term in this expansion corresponds to a sum of Feynman 
graphs that scale as a definite power of the expansion parameter i;<l. 

To leading order in velocity, the components of T M ^(x ,k) can be obtained from the 
diagram in Fig. HI The results are 

T 00 (x°,k) = ^m a e- 4kXa , (80) 
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(a) (b) (c) 

FIG. 5: Diagrams for v 2 corrections to T 00 . Dashed lines denote potential gravitons. 

which generates the moments 



rn 



J2m a + 0(v 2 ) (81) 

a 

X * = ^E m »< + (! V) ( 82 ) 
I» = jrm a [^i} TF + 0(v 2 ) (83) 

a 

l l3k = $> a [x* xjxj] TF + 0(^ 2 ) (84) 
/ d 3 xT° V [xV] TF = £ [xX] TF (85) 
where M = Y2 a m a- At order t> beyond the leading order we get the first contribution to 

- 1 - > 

T 0, (x°,k) = ^m a v:e- ik - x ", (86) 

which to this order in t> yields 

£>(v 3 ) (87) 

]m a (x a xv a ) l + £V) (88) 
;(v a x x a )X + (v a x x a )X) + 0(v 3 ) (89) 

/ rf 3 xT° V [xV] TF = ^ m a v a • x a [xX] TF (90) 

At order t> 2 beyond the leading order, there are corrections to T 00 from the graphs in 
Fig. as well as new contributions to T l i from the diagrams in Fig. For T 00 the results 
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(a) 



(6) 



FIG. 6: Diagrams needed to compute the leading contribution to T 4 - 7 . 
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(91) 
(92) 



2q 2 2(q + k) 2 m P/ 



(q 2 + q • k) e- ik Xb 



3 G N m a m h ^_. k . Xa 



a,b 



x a - x 6 



(93) 



with J q = J and where we have used the hHH Feynman rule derived in ref. [6]. To 
order v 2 we need the contributions of these graphs up to second order in the multipole 
expansion. This gives 



m 



a 

a 

a 

where we have defined a "renormalized" mass 



m a = m a 



1 2 1 V- 



G N m b 
x a - xJ 



(94) 
(95) 
(96) 

(97) 



To complete the full set of order v 2 corrections to the mass quadrupole P\ we also need 
to compute the leading order moments of T kk . We find 



Fig. E(a) = ^ m a v le~ 



«k-Xa 
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(98) 
(99) 
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and therefore 



a \ b ' a ' 



[<<] 



TF 



(100) 



Combining all of our results we finally obtain, from Eq. (!75|) . 



pj 



1 + 



I ^ d 



G N m b 



[«] TF + ^J2 m ^ ( x * i<<\ TF ) 



^ ma Jt ( Xa ' Va [ X « X ^ 



TF 



0( V 4 ). 
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It is straightforward to use the results obtained in this section to compute corrections to 
the quadrupole radiation formula. For illustration we consider two particles in a circular 
orbit of frequency Q, and work in the center of mass frame, X = 0, with X given in Eq. ( I95L 
Using the 1PN equations of motion, we obtain 

1 TF 



) x\ [xV] Ti? + ^x 2 (1 - 3i/) [vV]' 



•-/iv 7 ! - Av [(v x x)*x j + (v x x) J 'x*] 



(102) 

(103) 
(104) 



J*(f): 

r ifc (t) = //VI -4i/ [xVx*]' 

with x = xi — x 2 the relative coordinate, v = x, and \x = m^/M, v = fi/M. The PN 
expansion parameter is x = (GjvMO) 2 / 3 . The square modulus of the frequency domain 
moments (dropping terms involving S{oj) which do not contribute to radiation), 
vrT/iV r 



2Q 4 

2ft 4 
ttT^x 3 

aq 6 



i + 
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1 v I a; 

21 21 



[6(u - 2fi) + <?(w + 2fi)] , 



i - 4z/) - n) + 5(w + n)] 



(105) 
(106) 



l-4z/) 



5 



(107) 



with T = 27r5(0), can then be plugged into Eq. (fTB"|) . yielding the well known 1PN result 

po 



pO 



/1247 35 

1 — 1 V | X. 

V 336 12 



(108) 



Combining the results of this section with those of sec. IHH we can also obtain the 1.5PN 
and 2.5PN radiative (or "tail") corrections to the emitted power. This follows from Eq. ( IT21) 
and Eq. (I2"2"j) together with the 1PN mass monopole 

m v 



M 



1 



-x. 



(109) 
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One obtains 



Ptail * 3/2 f% 191 583 A 5/2 /1ln , 
-J22L = 47TX 3 / 2 _ 1 V 7TX H . (110 

Pgo V 672 24 / 

Here we have used the universality of the order rj 1 tail correction to obtain the coefficient 

of x 5//2 . Finally, including the results of sec. IIII B| we also obtain the term at 3PN that is 

non-analytic in the expansion parameter x 



Pi 



S56 x 3 lnx, (111) 
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where the renormalization scale has been chosen as //q ~ 1/r. 

VI. CONCLUSIONS 

In this paper, we have extended the EFT approach to compact gravitating systems, first 
proposed in p, to account for radiative corrections to gravitational wave emission. This 
includes corrections that arise from the multipole expansion a/ A <C 1, and from corrections 
to vacuum wave propagation, suppressed by powers of r\ = Gntti/X, that are present due 
to the non-linear nature of the Einstein equations. Many of the results obtained in this 
paper, in particular those of sec. IIII[ are universal, independent of the precise nature of the 
gravitating source. They can be used not only for post-Newtonian systems (as we consider 
in sec. |V]) but to other physical situations involving low frequency gravitons, such as soft 
bremsstrahlung in relativistic black hole collisions or graviton absorption by black holes. 



Of course, the 



approaches 



29 



FT approach necessarily shares some features in common with existing 
3l| to gravitational radiation from localized sources. For example, the sep- 

FT is physically equivalent to the standard matched 



aration of scales made explicit by an 

asymptotic expansions introduced in [28|, |29fl. However, there are some aspects of the prob- 
lem of gravitational radiation that are particularly well suited to EFT methods, namely 
the issue of how to handle the UV divergences that appear at sufficiently high order in 
perturbation theory (order rj 2 and beyond). Indeed, we saw that in the EFT, it is rela- 
tively straightforward to identify and isolate the terms in the perturbative expansion that 
contain UV poles, and to renormalize these terms by identifying a suitable counterterm in 
the effective theory Lagrangian. As we saw, this necessarily implies non-trivial RG flow 
of the multipole moments, which can be used to put constraints on the form of the non- 
analytic dependence on frequency of observables. Likewise, the identification, resummation, 
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and cancellation of IR divergences is made transparent in the way that the EFT computes 
perturbative corrections, as we discussed in sec. Ill II What emerges from our analysis is the 
following pattern of non-analytic terms in observables associated with radiation of frequency 
L>j from a fixed multipole 

oo 

\A e (u)\ 2 ~ S(u;)J2v 2n L [fi { £ + <D(r) 2 L) + 0( V 2 L) 2 + 0( V 2 L) 3 + •■■], (112) 

71=1 

where S(u) is the Sommerfeld factor in Eq. (I49j) . and L = In aw. The coefficient ft^n fixes 
the entire series of logs in the square brackets. 

There are several issues that have we not dealt with, but which we expect should be 
tractable by the methods presented here. For example, we have not included the effects of 
total angular momentum, which contribute beyond 3PN and were computed, using different 
methods, in 43]. It is clear that such terms are generated in our formalism by inclusion of 



the angular momentum coupling that appears in Eq. (CQ). We also have not considered in this 
paper the "memory" effect 43M45| which arises from the self interactions of the gravitational 
radiation. It corresponds to graphs with more than one coupling to multipole moments 
which can source radiation (quadrupole and higher) together with one or more graviton self- 
interactions. We do not see any obstruction to computing the resulting Feynman integrals by 
methods similar to those employed here. The matching to PN systems in sec. |V]has not been 
carried out to sufficiently high order to fix the matching scale /xo in the RG equation for the 
quadrupole moment. In principle, this computation is well defined in our approach, provided 
one regulates the theory consistently by dimensional regularization, as we have shown how 
to do here for the long distance i] corrections. Finally, we have restricted ourselves to 
time averaged observables generated by \A\ 2 . Thus we did not consider issues such as the 
computation of instantaneous waveforms or the back reaction of radiation on the multipole 
moments that serve as inputs into Eq. (CQ). It seems reasonable to imagine that some version 
of the formalism introduced in 13, [l4] could be applied to this problem. We hope to examine 
some of these questions in future work. 
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Appendix A: Feynman Integrals 

For the calculation of order rj 1 post-Minkowskian effects, the Feynman graphs in Fig. [1] 
can be expressed as integrals of the form (D = d — 1, with d the spacetime dimension) 

n - ln J (2^ (q 2)[( q + k ) 2 + (-k 2 -, e )]' 1 J 

with n < 4. Using rotational invariance, all such integrals can be converted ("tensor re- 
duced") to linear combinations of tensors constructed from k* and <5 lJ , with scalar integral 
coefficients of the form: 

d D q 1 



I (a) 



(2tt) c [ q 2]«[( q+k ) 2 + ( _ k 2_, e) ] 

1 r(a - D/2 + 1)T{D - 2 a -1) , 2 _ , D / 2 -a-! (A2) 



(4tt)d/2 T(D-a- 1) 



which has been performed by standard textbook methods, see for instance (46]. 

At order rj 2 , the loop graphs in Fig. |2]Ja),(b) are expressible as nested integrals, which 
can be done by first evaluating a subgraph proportional to the standard integral 



d D k 1 1 T{n l + n 2 -D/2)T{D/2-n 1 )T{D/2-n 2 ) 



p2\D/2-n 1 -n 2 



(27r)^[(k + p) 2 ] ni [kT 2 (4vr)^/ 2 r(m)r(n 2 ) T{D-m-n 2 

(A3) 

The result of this integration then yields linear combinations of terms of the form Eq. (1A1I) . 
which can be evaluated using Eq. (]A2[) . 

The graph in Fig. Wis) * s m ore challenging. Unlike Fig. |2]Ja),(b), it cannot be expressed 
in terms of simple nested one-loop integrals. After tensor reductions, all terms in Fig. Wis) 
can be expressed in terms of the scalar integrals (^=1,2,3 G 

/l({n i })(k) = / j—^ r/i ; x2ini rioin2 / /0 _: nr o1 n 3r , . ^2 . , ,0 — ( A4 ) 

l2(ni,n 3 ) = 



(2vr)^ [ (1 - k) 2 ] ni [P] " 2 J (2n)° [pa] "3 [ ( p + i) 2 + (_ k 2 _ ie) ] 

d D l 1 f d D p 1 



(2^ [ (1 - k ) 2 ] * [12 + (_ k 2 _ fe) ] J (2k)° [p2] [ (p + I) 2 + (_ k 2 _ fe) ] 

(A5) 
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The common integral over p can be calculated by standard Feynman parameter methods 
and expressed as a hypergeometric function. 

d D p 1 



hn(n 3 ) 



(^)° [p 2 ] n3 [(p + l) 2 + (-k2-*e)] 
1 T(n 3 - D/2 + l)r(£> /2 - n 3 ) ,_ k2 _ D /2-n 3 -x 



{Att) d / 2 r(D/2) 

l 2 



x 2 F 1 ln 3 - D/2 + 1, n 3 ; D/2; —j (A6) 

Alternatively (and equivalently, using a well known integral representation for the hyperge- 
ometric function), one can introduce a Mellin-Barnes representation 

- r * ni+zn _ t) j?+w , (A7) 



(p + l) 2 + (-k 2 -*e) 7- l0 o2« 7 (_ k 2_, e )-+i' 

where the integration contour runs parallel to the Imz axis in the way explained in [35]. In 
this representation 

hM = (4tt)^ r(n 3 ) ( " k 

+io ° rfz r(n 3 - d/2 + i + z)r(n 3 + ggx-g) / i 2 V 

_ 2™ r(D/2 + *) \-k*-ie) • 1 j 

Using Eq. ( 1A8j) together with Eq. ( 1A3j) . it is now possible to perform all momentum integrals 
in Eq. (1A4|) . with the result 

T fI X s 1 T(D/2-n 1 )T(D/2-n 3 ) , 2 . ^/2-n 3 -i /, 2 \0/2-ni-n 2 

'UWJ " (47r)D r(ni) fW l- k l k ) 



dz 

2~K% 



r(n 3 + z)T(l + n 3 - D/2 + z)r(£>/2 - n 2 + z) 

T(D/2 + z)T(D -n x -n 2 + z) 
r(ni +n 2 - z - D/2)T(-z] 



(A9) 



r(n 2 — 2) V ! + 

This integral can now be evaluated by residues. The result can be expressed in closed form 
as a sum over two generalized hypergeometric functions 4_F 3 evaluated at z — 1. It turns 
out that for the values of rij generated by the graph in Fig. 12(c), all such hypergeometric 
functions reduce to simple ratios of Gamma functions. As the expressions are best handled 
by computer algebra, we will not reproduce them here. 

For the second integral I 2 (rii,n 3 ), we use a different contour integral representation of 
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the hypergeometric function in Eq. (1A6|) (see Eq. D.72 of 351 ] ) 
1 1 



lining 



(A7r)D/2T(n 3 )T(D-l-n 3 : 
xT(D - 2n 3 - 1 - z)F(-z) 



dz 

too 27riL 
l 2 + (-k 2 - te 



r(l + n 3 - D/2 + z)T(n 3 + z) 



— k 2 — it 

Inserting this into Eq. (IA5[) . and performing the momentum integrals with the help of 

d D \ 1 



(A10) 



J (a,P) 



(2n) D [ 12 ]«[( 1 + k ) 2 + ( _ k 2_, e) ] 
1 r(a + P- D/2)T(D - 2a - 0) ,_^ 2 _ ^D/2-a-p 



(4tt)^/ 2 T{P)T{D - a - (3) 

we find the following contour integral representation for l 2 [ri2, %) 



(All) 



h{ri2,n 3 ) 



(4vr 



(-k 2 -,e) 



D— m— «3— 2 



X 



™ 27TZ 



r(n 3 )r(L»-l-n 3 ) 
r(l + n 3 - D/2 + z)r(n 3 + z)r(D - 2wi - 1 + z) 

T(D-n x -1 + z) 
T(D - 2n 3 - 1 - z)r(l + ni - D/2 - z)T(-z) 



(A12) 



r(i-*) 

For the case n 3 = 1, the contour integral can be done using Barnes' second lemma (see [35|, 
Eq. D.47), with the result 

1 (-k 2 - ie) D ^~ ni V{n x - D/2 + V 



D-3 (Aix) D V(D-m-l] 

r(m — D + 3)r(2D - 2m - 4)r(D/2 - ni)r(D/2 - 1) 



rp - 2m - i)r(2 - d/2) 

(A13) 



r(3D/2 - m - 3)r(m) 

For the case n 3 < 1, the contour integral is performed by taking residues. It turns out that 
only the pole at z = contributes, and the result is 



h(ni,n 3 < 1] 



_ k2 _ . u_ ni -ns-2 



(47T) 



x r(D - 2m - i)r(D - 2n 3 - i)r(m - D/2 + i)r(n 3 - d/2 + 1) 

x r(D - m - i)r(D - n 3 - 1) ' 1 j 

We did not bother to evaluate l2(n 2 ,n 3 > 1), as such terms are not needed to compute 
Fig. He). 
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Appendix B: UV poles 



After tensor reduction the amplitude for Fig. [2]^c) is a linear combination of scalar inte- 
grals of the form 

1 1 



(Bl) 



' 2 > " / 1>p [ (1 - k) 2 ] ni [P] - [l 2 + (-k 2 - ie) ] [ P 2] - [ (p + l) 2 + (-k 2 - ie) ] 

with m, ri2, n 3 G Z. Here we have defined f p = f d D p/(2ir) D . As in appendix lAl D = d — 1. 
In order to extract the UV behavior of the integrals we formally expand the denominator in 
powers of — k 2 — ie, which yields 

oo 

J 123 = (_i)^'(_k 2 - te) l+j / --n • (B2) 

i,i= ■ /i .p L l 1 - k ) J L 1 J Lp J L (p + ! ) J 

The p-integration is then easily performed using Eq. ( ]A3j) 

r Vfl:^+» 1 r(n 3 +j + l-D/2)r(D/2-n 3 )r(D/2-j-l) 

123 J (47r)f r(n 3 )r(j + i) r p-n 3 -i-i) 

x f w, (B3) 

Ji [(i + k) 2 ] ni \p~y> 2+n ' i+i+j+2 ~~ 

and we note that the first line of Eq. ( TB3I) is free of any divergences for d = 4. The remaining 
task is therefore to extract the UV poles from the integral over 1 in Eq. ( 1B3I) which is 8 

UV divergent if n\ + n 2 + n 3 + i + j + 3 < d, (B4) 
IR divergent if n 2 + n 3 + i + j + 3 > d . (B5) 

Only a finite number of combinations of i and j can yield UV divergences which truncates 
the sums in Eq. ( 1B3|) . If n\ < it is possible that the integral is both IR and UV divergent, 
and our computation involves the cases n\ = 1, 0, —1, —2, —3. 

In the case of n\ — 1 the UV divergences never overlap with the IR divergences and we can 
extract the UV divergences using Eq. ( 1A3I) and keeping only the UV poles. For n% > the 
integration in Eq. flB3j) is a scaleless integral which vanishes in dimensional regularization. 
When it is logarithmically divergent however the vanishing implies a cancelation of IR and 



Here we took into account that ni < 1 for the integrals needed so that the first factor in the denominator 
never yields an IR divergence. 
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UV divergences. In order to extract the UV poles, we only need to consider this case and 
use 



1 



1 



' 47T 2 (d-4V 

In this way we find for the UV divergent part of diagram Fig. [2{c) 



A {c) 



(G N m\k\Y 



uv 



(k 2 + te) 



7r/r 



(d-4) 



1046 1 
315 (d-4) 



(B6) 



(B7) 
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